Prestack first-break traveltime tomography using the double-square-root eikonal equation

Traveltime tomography with shot-based eikonal equation fixes shot positions then relies on inversion to resolve any contradicting information between independent shots and achieve a possible cost-function minimum. On the other hand, the double-square-root (DSR) eikonal equation that describes the whole survey, while providing the same first-arrival traveltimes, allows not only the receivers but also the shots to change position and therefore leads to faster convergence in tomographic inversion. The DSR eikonal equation can be solved by a version of the fast-marching method (FMM) with special treatment for its singularity at horizontally traveling waves. For inversion, we use an upwind finite-difference scheme and the adjoint-state method to avoid explicit calculation of Fréchet derivatives. The proposed method generalizes to the 3D case straightforwardly.